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William Gelbach was the Director. Mr. James W. Adams was the 
technical manager. Most of the analytical effort in support 
of this project was provided by Dr. Frederick Chichester*, 
who wrote all sections of this report, and Mr. Don Lipski of 
the Bendix Simulation Center. Computer simulation support 
was provided by Mr. Isaac Emmanuel, Mr. Don Lipski and Mr. 
Alexander Labounsky of the Simulation Center. The guidance 
of Dr. Henry B. Waites of MSFC during the course of this 
study is gratefully acknowledged. 


♦Presently a member of the Technical staff of the Guidance 
Systems Division, Teterboro, New Jersey. 


ABSTRACT 


Application of modular control techniques to the attitude 
control of a prototype flexible spacecraft and a prototype 
flexible space platform was further developed by determining 
numerical values for the physical parameters of a four body 
approximation of the MSFC/hybrid deployable truss incorporat- 
ed in the space platform model, generating sensitivity coef- 
ficients for the model of the flexible spacecraft, evaluating 
the changes in the digital computer simulation of the flexi- 
ble spacecraft resulting from the addition of another rigid 
body to the model and comparing attitude control effective- 
ness with actuators on more than one rigid body of the model 
with that for the case in which the actuators were restricted 
to one body. 
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SECTION 1 


1.0 INTRODUCTION 

This report is submitted in compliance with the Scope of Work 
under contract NAS8-33979. The period of performance covered 
by the contract is from October 15, 1981 to October 31, 

1982. The submission and approval of this report constitute 
the successful completion of the "Exhibit B" portion of the 
contract. 

This report is a sequel to three others, two of them previ- 
ously submitted under a different contract number. The two 
prior reports, under a different contract number, references 
(1-1) and (1-2), were submitted in October 1978 and Septem- 
ber, 1979 and covered the periods from July 27, 1977 to July 
27, 1978 and from August 26, 1978 to August 26, 1979, respec- 
tively, in compliance with "Exhibit A" of contract NAS 8- 
32660. 

The prior report under contract NAS-33979, reference (1-3), 
was submitted on March 8, 1982 and covered the period from 
August 15, 1980 to October 15, 1981 in compliance with "Ex- 
hibit A" of the contract. 

1 . 1 OBJECTIVES 

The sections that follow summarize the effort expended on the 
Modular Design Attitude Control System Study contract. The 
overall objective of the study was to further develop a new 
approach to applying attitude control to mathematical models 
of the rotational dynamics of a representative flexible 
spacecraft and of a representative space platform incorporat- 
ing both the flexible spacecraft and a deployable truss. 
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More specifically, to complete the digital computer simula- 
tion of the ten body model of the prototype flexible space 
platform initiated during the previous reporting period, it 
was necessary to determine numerical values for the physical 
parameters of the four body approximation of the MSFC/hybrid 
deployable truss that had been developed algebraically. Ad- 
ditional specific objectives included development of numeri- 
cal values for the relative sensitivity of the state varia- 
bles of the flexible spacecraft model to changes in the nu- 
merical values of its physical parameters, evaluation of 
changes in the digital computer simulation of the model re- 
quired to accommodate addition of another rigid body and com- 
parison of attitude control effectiveness between the case in 
which actuators were distributed over more than one rigid 
body of the model and the case in which they were restricted 
to one body. 

1.2 SCOPE 

Study effort was concentrated in five main areas: 

1 . Development of a dynamically equivalent four body approx- 
imation of the NASTRAN finite element model supplied for 
the MSFC/hybrid deployable truss to support the digital 
computer simulation of the ten body model of the flexible 
space platform that incorporates the four body truss mod- 
el . 

2. Generation of coefficients for sensitivity of state vari- 
ables of the linearized model of the three axes rotation- 
al dynamics of the prototype flexible spacecraft with 
respect to the model's parameters. 
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3 . 


Evaluation of software changes required to accommodate 
addition of another rigid body to the five body model of 
the rotational dynamics of the prototype flexible space- 
craft. 

4. Comparison of effectiveness of attitude control for actu- 
ators on two bodies of the six body model of the proto- 
type flexible space platform with that for actuators re- 
stricted to one body of the same model. 

1 . 3 GENERAL 


This report is comprised of six sections. Section 2 de- 
scribes the determination of the numerical values of the pa- 
rameters of the four body model of the MSFC/hybrid deployable 
truss in such a way that five modes of its rotational dynam- 
ics closely matched those of NASTRAN finite element model de- 
scribed in Ivey (1-4). Section 3 describes the generation of 
coefficients of the sensitivity of the state variables to 
changes in the parameters of the linearized five body model 
of the rotational dynamics of the prototype flexible space- 
craft. Section 4 presents an evaluation of the changes in 
the computer simulation of the five body model of the proto- 
type flexible spacecraft required to accommodate the addition 
of another rigid body to the model. Section 5 presents a 
comparison of the effectiveness of modular attitude control 
for actuators on two bodies of the six body model of the pro- 
totype space platform with that for actuators restricted to 
one body of the same model. Section 6 lists a number of con- 
clusions and recommendations drawn from the results of the 
tasks described above. References are listed at the end of 
each section. 
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The original RFQ requested that the International System of 
units (designated as SI) be used in the program and in any 
reporting. Torques, moments, angular momentum, moments of 
inertia and distances, however, are stated in English units 
since this was the system of units used in presenting all of 
the vehicle data in the RFQ and the truss data in Ivey (1-4). 

1.4 REFERENCES 

1-1 Guidance Systems Division, The Bendix Corporation, "Space 
Construction Base Control System", Final Report, Contract 
NAS8-32660 for George C. Marshall Space Flight Center, 
October 27, 1978. 

1-2 Guidance Systems Division, The Bendix Corporation, "Space 
Construction Base Control System", Final Report, Contract 
NAS8-32660 for George C. Marshall Space Flight Center, 
September 1, 1979. 

1-3 Guidance Systems Division, The Bendix Corporation, "Modu- 
lar Design Attitude Control System", Final Report, Con- 
tract NAS8-33979 for George C. Marshall Space Flight Cen- 
ter, March 8, 1982. 

1-4 Ivey, Wayne, "Vibration Analysis of the MSFC/Hybrid De- 
ployable Truss", Marshall Space Flight Center memorandum, 
June 26, 1980. 
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SECTION 2 


2.0 APPROXIMATING A NASTRAN FINITE ELEMENT MODEL OF 
THE HYBRID DEPLOYABLE TRUSS BY A FOUR BODY MODEL 

2 . 1 INTRODUCTION 

In Chichester (2-1), a prototype flexible spacecraft com- 
prised of a rigid central body with two flexible solar wings 
was approximated by an assembly of five rigid bodies inter- 
connected by a spring hinge suspension shown in Figure 2-1 . 
The resulting linearized model of the rotational dynamics of 
the spacecraft was written in state variable form to facili- 
tate the application of multilevel attitude control. In lat- 
er work, Chichester (2-2), this prototype flexible spacecraft 
was connected to another spacecraft via a hybrid deployable 
truss to form a flexible space platform. This flexible space 
platform was approximated by the assembly of ten rigid bodies 
interconnected by a spring hinge suspension, depicted in Fig- 
ure 2-2, under the assumptions that the second spacecraft 
could be approximated by a single rigid body and the truss 
could be approximated by four rigid bodies serially connected 
by a spring hinge suspension. The rotational dynamics of the 
MSFC/hybrid deployable truss were presented in Ivey (2-3) in 
terms of a NASTRAN finite element analysis of the response of 
the truss to perturbations while it was in a cantilever con- 
figuration. 

The work reported in the present section was directed toward 
approximating the NASTRAN model of the truss with a model 
consisting of four bodies serially connected by a spring 
hinge suspension. This approximation was effected by appli- 
cation of the following procedure: 
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FIGURE 2-1. 

TOPOLOGICAL DIAGRAM OF FIVE BODY 
APPROXIMATION OF A PROTOTYPE FLEXIBLE SPACECRAFT 
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FIGURE 2-2. 

TOPOLOGICAL DIAGRAM OF TEN BODY APPROXIMATION 
OF TWO SPACECRAFT INTERCONNECTED BY A DEPLOYABLE TRUSS 
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1 . Decompose the truss into four serially connected modules 
determining the total mass and rotational inertias for a 
rigid body with its mass center on the axis of the truss 
for each module. 

2. Develop five body cantilever truss model. 

3. Develop undamped form of cantilever truss model. 

4. Generate spring coefficient values to match eigenvalues 
of undamped five body cantilever truss model with those 
of the NASTRAN model. 

5. Generate damping coefficients to add n% of critical damp- 
ing to the five body model of the cantilever truss with 
eigenvalues matching those of the NASTRAN model. 

2.2 DECOMPOSITION OF TRUSS INTO MODULES 

A perspective drawing of the MSFC/hybrid deployable truss 
that appeared in Ivey (2-3) is reproduced in Figure 2-3 with 
shaded circles representing the locations of the centers of 
mass of four 8000 pound payloads rigidly attached to the 
fifth and tenth cubic submodules of the truss. The truss and 
its payloads were decomposed into two pairs of identical mod- 
ules serially connected by three degree of freedom hinges 
shown in Figure 2-4. Each module in one pair consists of 
only four cubic submodules of the truss. Each module in the 
other pair consists of one cubic submodule of the truss with 
8000 pound payloads rigidly attached both above and below 
it. An equivalent rigid body with its center of mass located 
on the centerline of the truss and the same total mass and 
rotational inertias was determined for each of the four mod- 
ules under the following assumptions: 
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FIGURE 2-3. PERSPECTIVE VIEW OF HYBRID DEPLOYABLE TRUSS 
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FIGURE 2-4 

SIDE VIEW OF MODULES OF DECOMPOSED TRUSS 
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1 . The mass of each cubic truss submodule is uniformly dis- 
tributed m the cube. 

2. The mass of each payload may be represented as a point 
mass at its center of mass. 

Details of the analysis to determine these equivalencies ap- 
pear in Chichester and Emmanuel (2-4). 

The equivalent total masses and rotational inertias for the 
rigid bodies representing the modules are presented in Table 
2-1. The parameters utilized in this table are the follow- 
ing : 

m 1 = total mass of cubic submodule of truss 
m 2 = total mass of 8000 pound payload 

L 1 = 1/2 length of one side of cubic submodule of truss 

L ? = distance of center of mass of payload from nearest 
face of cubic submodule of truss. 

The determination of the total masses and rotational inertias 
of the equivalent rigid bodies reduced the model of the truss 
to four rigid bodies with their centers of mass located on 
the axis of symmetry of the truss serially connected by three 
degree of freedom hinges with associated rotational spring 
and damping coefficients. 

.3 DEVELOPMENT OF CANTILEVER TRUSS MODEL 

The vibrational modes obtained for the finite element NASTRAN 
model of the hybrid truss in Ivey (2-3) resulted from per- 
turbing the end of the truss to which the payloads are 
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TABLE 2-1 . 

EQUIVALENT TOTAL MASSES AND 
ROTATIONAL INERTIAS OF TRUSS MODULES 



MODULES W/O 
PAYLOAD 

MODULES WITH 
PAYLOAD 

Total Mass 

4n^ 

2ni2 + 

J x 

8 . 2 
3°1 L 1 

n 2 [ (L x + L 2 ) 2 + h 2 ] 

*y 

68 T 2 

sane as I x 


sane as Iy 

2™ L 2 

3°1 L 1 
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attached while its unloaded end was attached to a wall in a 
cantilever configuration. To replicate these conditions, the 
unloaded end of the four body model of the truss was attached 
via a three degree of freedom hinge to a fifth body of large 
mass and rotational inertia with its center of mass close to 
this hinge. The imposition of these conditions yields re- 
sponses to displacements of the free end of the truss that 
very closely approximate cantilever responses according to 
investigations by Lipski (2-5). 

The topological diagram of the five body cantilever truss 
model is depicted in Figure 2-5 which is of the same form as 
the one developed for the prototype flexible spacecraft in 
Chichester (2-1) but with different values for the masses and 
rotational inertias of the rigid bodies. From equation (19) 
of Chichester (2-1) the linearized state variable representa- 
tion of the five body rotational dynamics model of the truss 
without either external or control torques may be written as 
follows : 

. 

a) = GLC a) + GLK ct 

— S— s— 

A 

a = K(£ 
where : 

T T T.T 

(a) ( (i) j U) / « • • f (0 ) 

“ “1 2 5 

w iy , “> iz ) = rate of angular rotation of ith 
rigid body 

"T "T "T.T 

a =(a, a,..., a) 

- _ 2 _ 5 


a = a a = a a = a a = a a = a 

— 1 — 2 4 ~2 —3 5 —3 —12 ~ 4 “13 ~ 5 “ 1 


( 2-1 ) 
( 2 - 2 ) 
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FIGURE 2-5. 

FIVE BODY MODEL OF TRUSS IN CANTILEVER CONFIGURATION 
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a = (<p , 0 , i|> ) 
— 1 ill 


a . .= (<)>.., 0 . . , ij> . . ) 
— l] i] l] 1] 


<1> l / 9 x r ^ 


= inertially referenced Euler angles of 
body 1 


't’l]' e ij' ^i] 


= relative Euler angles of body j with 
respect to body i 


G, L, C s , K s , K = coefficient matrices defined in 

Chichester ( 2-1 ) 


2.4 DEVELOPMENT OF UNDAMPED FORM OF THE CANTILEVER TRUSS MODEL 


With the mass and rotational inertias determined for each 
rigid body in the five body model of the cantilever truss the 
values of the spring and damping coefficients associated with 
the three degree of freedom hinges were determined in such a 
way that the first several vibrational modes of the five body 
model closely matched those of the NASTRAN finite element 
model. The basis for this matching was the adjustment of 
these parameters until the corresponding eigenvalues of the 
two models were closely matched. To this end, the five body 
cantilever truss model was transformed to its undamped 
form. This was accomplished as follows. Differentiating 
equation (2-2) with respect to time yields: 

a = Ku (2-3) 

Premultiplying both sides of equation (2-2) by K -1 yields: 

a 

w = K~ (2-4) 
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Substitution of equations (2-3) and (2-4) into (2-1) yields 
the basic equation for the unforced system in terms of Euler 
angles . 

* "a 

a = KGLC K -1 a + KGLK a (2-5) 

— s — s 

The undamped form of the unforced system model was written by 
omitting the term including 

•• <*> __ 
a - KGLK a = 0 (2-6) 

— s — — 

The undamped natural frequencies corresponding to the modes 
of the system described by equation (2-6) are the square 
roots of the corresponding eigenvalues of the system as shown 
m Lipski (2-5). 

The above undamped model was partitioned according to spatial 
axis as follows: 

± ~ (KGLK s ) x ± = 0 

0 - (KGLK ) 0 = 0 (2-7) 

— ' s y— — 

\Ji - (KGLK ) = _0 

s z 

where : 


1 = 

( <t> f 

♦ / 

<t> , 

<f> r 

♦ ,) 


24 

3 b 

12 

1 3 

1 

0_ = 

(9 , 

9 F 

0 , 

9 

0 ) 


24 

3 5 

12 

13 

l 

1 = 

( V r 

'I' f 

'f’ f 

'J' , F 

* ) 


24 

3 5 

12 

1 3 

l 
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KGLK = 
S 


( KGL K s ) x 

[ 0 ] 

[ 0 ] 


[ 0 ] 

(KGLK s ) 

[ 0 ] 


[ 0 ] 

[ 0 ] 

(KGLK ) 

' s z 


£ = null vector of dimension 5 

[0] =5x5 null matrix 

2.5 GENERATION OF SPRING COEFFICIENTS TO 
MATCH EIGENVALUES OF UNDAMPED MODELS 


The undamped five body cantilever truss model was simulated 
on a digital computer to generate values of the spring coef- 
ficients that correspond to close matching between eigenval- 
ues of the five body model and those of the NASTRAN finite 
element model. Since the five body cantilever truss model 
was based upon decomposition of the truss itself into four 
modules, matching was sought for just five of the eigenvalues 
corresponding to two vibrational modes about each of two of 
the spatial axes and one about the third axis. To attain 
this matching an algorithm utilizing a steepest descent or 
gradient method was programmed to systematically vary the 
values of the 12 spring coefficients until the five eigenval- 
ues of the five body model of the undamped cantilever truss 
were closely matched with those of the NASTRAN model. 

Implementation of this approach resulted in matches between 
corresponding eigenvalues of the two truss models with an 
accuracy of three decimal digits. The corresponding values 
of the spring coefficients are displayed in Table 2-2. 
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TABLE 2-2. 

SPRING COEFFICIENTS 

CORRESPONDING TO HATCHING OF FIVE EIGENVALUES 
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2.6 


DETERMINATION OF DAMPING COEFFICIENTS 
CORRESPONDING TO n% OF CRITICAL DAMPING 


The vector matrix equation representing the damped unforced 
five body model of the cantilever truss, equation (2-5), was 
partitioned according to spatial axis by the same procedure 
as was utilized for the spatial partioning of the undamped 
model. The resulting set of three vector matrix equations 
was then expanded to a set of 12 scalar equations by writing 
the first four scalar components of each vector matrix equa- 
tion corresponding to a spatial axis. Each of these scalar 
equations corresponded to one of the relative Euler angles, 

4> -j for the x axis, 0^ for the y axis, and for the z 

axis, where i] = 24, 35, 12, 13. In general, the equation 
for each relative Euler angle contained terras representing 
coupling with the other three relative Euler angles about the 
same axis. The coupling terms were ignored and the Laplace 
transform was applied to the resulting set of decoupled sca- 
lar equations. Since, m general, (s), 0 x ^(s) and 
¥ x -j ( s ) * 0, their coefficients in the frequency domain must 
then be equal to zero resulting in the following set of 24 
equations . 


xaxis: s 2 + b . s+c =0 


1]X 


13* 


b c 

1 px _ SI] x 

c k . 

1JX S1]X 


y axis: s 2 +b s+c =0 

13Y iDY 


b c 

A3Y = 

c i]y k si]y 


( 2 - 8 ) 


z axis: s 2 +b s+c- =0 


1JZ 


13 z 


b . c 

ljz _ sij z 

c . k 

1JZ sijz 


where : 


lj = 24, 35, 12, 13 
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For critical damping, the discriminant of each of the above 
quadratic equations equals zero. Applying this condition led 
to the following set of equations for the damping coeffi- 
cients in terms of the spring coefficients corresponding to 
n% of critical damping. (Damping coefficient values obtained 
for n = 1 appear in Table 2-3.) 


c 

S1]X 

c 

sijy 

c 

SI] z 


• 02n(o ijx ) 


.02n(o ijy ) 


•02n(c . ) 

' 1] z 


J /2 

_1/ 2 


si jx 
k sijy 

^Sl j z 


(2-9) 


Details on the matching of the rotational dynamics of the 
four body model to those of the NASTRAN model of the truss 
appear in Chichester and Emmanual (2-6). 
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TABLE 2-3. 

DAMPING COEFFICIENTS CORRESPONDING 
TO 1% OF CRITICAL DAMPING WITH MATCHED EIGENVALUES 


£ ft. -lb /rad /sec.J 
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SECTION 3 


3.0 GENERATING NUMERICAL VALUES FOR SENSITIVITY COEFFICIENTS OF 
SCALAR STATE VARIABLES WITH RESPECT TO MODEL PARAMETERS 

3 . 1 INTRODUCTION 


Earlier work in the application of modular control to flexi- 
ble spacecraft, Chichester (3-1), (3-2) and (3-3), was predi- 
cated upon the assumptions that all of the numerical values 
of the physical parameters of the model of the flexible 
spacecraft being controlled were known to sufficient accuracy 
and were sufficiently time-invariant to support effective 
attitude control. Since the physical parameters of actual 
spacecraft rarely satisfy these rather broad assumptions, it 
is necessary to assess the sensitivity of the time trajecto- 
ries of the state variables of the rotational model of the 
flexible spacecraft to small changes in the numerical values 
of each of these parameters. 

The study presented in this section is based upon the linear- 
ized state variable model of the rotational dynamics of the 
five body approximation of the prototype-flexible spacecraft 
developed in Chichester (3-4) and depicted in Figure 2-1. 

The optimal attitude control problem associated with this 
model may be stated as follows: 


(_x T (3x + u^W^Jdt (3-1) 


Minimize : 


- / 1 


3-1 


subject to the linearized state variable (constraint) equa- 
tion: 

x = Ax + Bu; te[t Q , t f l (3-2) 


where : 

x = (u> T , a^)^ = state vector of dimension 30 

«0 = (u> T , u/ r , U) T ) T = body angular rate vector of 

- — 1 —2 — 5 . . , n 

dimension 15 


T.T 


~ T T T T T x 

° = ( — 2 4 ' -05' -12' -13' -1 } 


co k = (ai kx , » kyf u> kz ) T = angular rotation rate of kth body 


-ij (<f> ijx' 9 ijy' ^ijz^ 

'f’ij' 9 ij' '•'ij = relative Euler angles of body j with 

respect to body i 


<j> i r e tj; i = inertially referenced Euler angles of 
central body (body 1) 

Q = state vector error coefficient matrix of dimension 
30 x 30 


W = control vector energy coefficient matrix of 
dimension r x r (r _<_ 15) 

A = state vector coefficient matrix of dimension 
30 x 30 


3-2 



B = control vector coefficient matrix of dimension 
30 x r 

= control vector of dimension r 
t Q = initial time 
tf - final time 
3 . 2 PROCEDURE FOLLOWED 

1 . Express the state vector coefficient matrix, A, and the 
control vector coefficient matrix, B, in the state varia- 
ble linearized model of the rotational dynamics of the 
five body approximation of the prototype flexible space- 
craft in terms of the model parameters p j , j = 1, 2, ..., 
54 which represent: 

o five body masses 

o 15 non-zero coefficients of rotational inertia 

o 10 non-zero spatial components of location vectors for 
mass centers of the rigid bodies comprising the model 

o 12 spring coefficients 

o 12 damping coefficients 

2. Corresponding to the perturbed value, 

P'j = P D + Apj, (3-3) 

of each model parameter develop a perturbed state varia- 
ble model of the rotational dynamics of the flexible 
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spacecraft by incorporating the corresponding perturbed 
coefficient matrices. 


= A 


p . + p . 


and B" = B 


P- Pj 


(3-4) 


3. Develop the corresponding perturbed optimal attitude con- 
trol problem for each perturbed model generated in Step 
2. For each model parameter perturbation, Apj, the cor- 
responding perturbed optimal attitude control problem is 
expressed as follows: 

Minimize : 


P' = / tf l(x') T Qx' + (u') T W u u']dt (3 _ 5) 

fc o 

subject to the perturbed linearized state variable (con- 
straint) equation, 

x' = h'x' + B'u'. (3-6) 

The solution vectors for this control problem 
are _x'(t) and u'(t) for te[t Q , t^], 

4. Simulate each perturbed optimal attitude control problem 
of Step 3 on a digital computer and obtain the perturbed 
control vector, response u'(t), to a standard initial 
displacement for each perturbed model parameter, P'j. 

5. Utilize a digital computer simulation of the unperturbed 
state variable model (with coefficient matrices A and B 
instead of A^ and B' , respectively) to generate the time 


3-4 


trajectories of the state vectors, x.(t) • °f the unper- 
turbed model to each of the perturbed control vector time 
trajectories, u^(t), generated in Step 4 in conjunction 
with the same initial displacements. The unperturbed 
linearized state variable model with the perturbed con- 
trol vector, u^(t), as a forcing function then takes the 
following form 

x" = Ax" + Bu' (3-7) 


where x" (t) represents the resulting state vector solu- 
tion for te[t 0 , tf]. 


6. Use scalar state responses, x^(t), of the original model 
without any perturbations in conjunction with the corre- 
sponding scalar components of the state vectors obtained 
in Step 5 to generate the following scalar state pertur- 
bations . 


A Xi (t) = xt'(t) - x.(t) te[t Q , t f j 


(3-8) 


7. Construct a table of maximum magnitude sensitivity coef- 
ficients, S i j m , by utilizing the following definition. 


Max 

S ijm = te ^o' 


AXi 

f J 


(3-9) 


8. Construct a table of "steady state" sensitivity coeffi- 
cients, S^jgg, by using the following relationship 

Ax i(t f ) 


S ijss 


Ap, 


(3-10) 
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3 . 3 RESULTS 


Representative numerical results of the application of the 
procedure described above are displayed in Tables 3-1 and 3-2 
with the entry of largest magnitude in the upper left corner 
and the remaining entries in descending order across the rows 
and down the columns. The entries of Table 3-1, the maximum 
magnitude sensitivity coefficients, S^j m , provide an index to 
the relative sensitivity of the scalar state variables, x^, 
to changes in the numerical values of the model's parameters, 
p j , while the controlled system is responding to a distur- 
bance. The entries of Table 3-2, the steady state or final 
time sensitivity coefficients, S^j ss , provide an ordering of 
the relative sensitivity of the state variables, x^, to 
changes in the numerical values of the model's parameters, 

Pj, when the controlled system is at an equilibrium condi- 
tion. In both tables the Euler angles 4> , 9, \J>, represent 
rotations about the x, y and z spatial axes, respectively, 
and the double subscript, i j , on the relative Euler angles 
represents angular displacement of body j with respect to 
body i. 


Comparison of the two tables reveals certain patterns common 
to the two. For example, the six largest transient sensivity 
ratios are, in descending order along the first column of 
entries m Table 3-1 , 


A *i | 

A c * 
sj bZ max 


A<j> L j 

Av i*z max 


(3-11) 


* v ijz represents the z axis component of the location vector 
of the jth hinge with respect to the mass center of the ith 
body. 
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The same ordering applies to the six largest steady state 
sensitivity ratios in the first column of entries of Table 
3 - 2 . 


Although the detailed ordering of the sensitivity coeffi- 
cients differs between the two tables beyond this point, cer- 
tain common patterns may be discerned. The changes in model 
parameters to which the state variables are most sensitive 
are, in descending order, changes in the z axis damping coef- 
ficients, c SJ5z , c S24Z , c sl3Z , c sl2Z , and changes in the z 
axis components of the location vectors of the hinges con- 
necting body 1 to the adjacent bodies, v liz and v 12Z . In 
addition, among the seven scalar state variables showing the 
highest sensitivity to model parameter changes in both ta- 
bles, six, the inertially referenced angular displacements of 
body 1, $ , 6 L and \|> lf and the relative angular displacements 
between the bodies close to the free ends of the appendages, 
i)> 24 , and $ ib , are common to both although they appear 

with different ordering. 
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SECTION 4 


4.0 EVALUATING SOFTWARE CHANGES REQUIRED TO ACCOMMODATE 
ADDITION OF ANOTHER RIGID BODY TO AN EXISTING DISCRETE 
MASS MODEL OF A FLEXIBLE SPACECRAFT UNDER MODULAR 
ATTITUDE CONTROL. 

4.1 INTRODUCTION 


Earlier analyses of linearized state variable models of the 
rotational dynamics of flexible spacecraft under modular at- 
titude control were predicated upon models comprised of a 
constant number of flexibly interconnected rigid bodies as in 
Chichester (4-1), 4-2), (4-3) and (4-4). However, situations 
can be anticipated in which it would be necessary to add one 
or more rigid bodies of appreciable weight to a flexible 
spacecraft while it is in orbit. This could occur, for exam- 
ple, in docking one spacecraft to another. For this reason, 
it was deemed desirable to evaluate systematically the chang- 
es in a representative flexible spacecraft model that would 
result from adding another rigid body to it. Furthermore, an 
attempt was to be made to recast the model of the flexible 
spacecraft into a form that would experience minimal changes 
due to the addition of a rigid body. 

4.2 PARTITIONING OF ATTITUDE CONTROL 
PROBLEM FOR FIVE BODY MODEL 


1 . The linearized state variable model of rotational dynam- 
ics of the three axis five body approximation of the pro- 
totype flexible spacecraft depicted in Figure 2-1 was 
partitioned according to spatial axis and the y axis sub- 
model was found to be decoupled from the submodel for the 
x and z axes. 
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2. The state vectors of each of the submodels of Step 1 was 
partitioned into a rigid body angular rate vector and an 
Euler angle vector for each spatial axis. 


3. Each of the three spatial body angular rate vectors and 
the three spatial Euler angle vectors of Step 2 was par- 
titioned into components corresponding to the five rigid 
bodies comprising the model. The resulting pair of par- 
titioned submodels was written in the following form. 


xz axes submodel: 


x = A x + B u 

— l i—i i—i 


(4-1) 


y axis submodel: 

x = A x + B u 

—i 2—i 2~2 


(4-2) 


where : 


x = (iu r , 4 >' r , (o T , 4> T ) T = xz axes state vector 
I X z 

T T T 

x = (^y, £ ) = y axis state vector 

m T T 

u = (u , u ) = xz axes control vector 

— 1 —X — z 


u = u = y axis control vector 
-x -y 


(4-3) 

(4-4) 

(4-5) 

(4-6) 


— X 

= 

(U 1X' 

U 2X' 

• 

• • 

r 

W 5X 

^z 

= 

< u iz' 


• 

• • 

9 

“5Z 

3| 

= 

(“iy ' 

>i 

N 

3 

• 

• • 

t 

“sy 

i 

= 

(♦ r 

4> f 



f 

<P 



24 

3 b 


12 


1 3 

i 

= 

('I' r 

V r 

'P 


r 

* 



24 

3 5 


1 2 


1 3 

8 

= 

(9 , 

9 , 

9 


r 

0 



24 

3b 


12 


1 i 


) T 

) T 

) T 


f 4» ) 

i , 

/ 'P , ) 

1 r 

, 9 ) 

1 


(4-7) 


(4-8) 
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“x = 

< U iX' 

U 2X' 

• • • f U rx ) } 

r 

1 

Mz = 

< U iZ ' 

U 2Z' 

11 \ T. 

• • • / u rz ' 9 

m 

r 


|C 

►< 

II 

(U iy , 

u zy' 

• • • t u r y) 9 

r 

<. 


= xz axes state vector coefficient matrix of dimen- 
sions 20 x 20 for 5 body model. 

= y axis state vector coefficient matrix of dimensions 
10 x 10 for 5 body model. 

B = xz axes control vector coefficient matrix of dimen- 
sions 20 x r for 5 body model. 

B = y axis control vector coefficient matrix of dimen- 
sions 10 x r for 5 body model. 

4. The optimal attitude control subproblem associated with 
jth submodel (j = 1, 2) was written as follows: 

Minimize : 
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subject to: 


- Vj + Vd 


where: 


= jth state vector 


x d = jth desired state vector 


(4-11 ) 
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= jth control vector 


Qj = jth state vector error weighting coefficient matrix 
Wj u = jth control energy weighting coefficient matrix 
t Q = initial time 
tf = final time 


5. Application of linear quadratic regulator (LQR) tech- 
niques presented in Chichester (4-5) reduced each of the 
optimal attitude control subproblems to the solution of 
Riccati equations of the following form. 


K. 

D 


-(K. A. 
1 3 


A T K 


K.RK. 
3 3 


V 


j - 1 


(4-12) 


where : 


Rj = -B j^j u~ 1 Bj'*’ (4-13) 

K j ( t f ) = [01 (4-14) 

For j = 1 , the 5 body xz axes subproblem Kj is a 20 x 20 
symmetric matrix? for j = 2, the 5 body y axis subproblem 
is a 10 x 10 symmetric matrix. 

4.3 PARTITIONING OP THE SIX BODY MODEL ATTITUDE CONTROL PROBLEM 


The preceding steps were applied to the linearized state va- 
riable model of the rotational dynamics of the six body con- 
figuration depicted in Figure 4-1. 

The state variable model for the six body configuration was 
written with the same structural form as that shown for the 
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SPACECRAFT 

1 



FIGURE 4-1 

TOPOLOGICAL DIAGRAM OF SIX BODY MODEL 
OF TWO SPACECRAFT INTERCONNECTED BY A DEPLOYABLE TRUSS 
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five body configuration in equations (4-1) and (4-2) with the 
connection of body 6 to body 1 being accommodated by incre- 
menting the dimension of the state subvectors, u» , w , co , 
and 9 f and the maximum allowable dimension/ r f of the con- 


. vectors u 

x' 

and 

u , by one 

-y 

as 

follows. 


0) 

— X 

= 

(w ix» 

“zx' 

• • • f 

“sx' td bX ) 
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• • • f 
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(4-17) 

% 

= 

( U iy f 

u zy ' 

• • « f 

u ry ) T ; r 

1 6 




The rigid body angular rate subvectors of equation set (4-15) 
and the Euler angle subvectors of equation set (4-16) were 
aggregated into an xz axes submodel state vector and a y axis 
submodel state vector for the six body model of the same form 
as that displayed for the five body model in equations (4-3) 
and (4-4) and the control subvectors of equation set (4-17) 
were aggregated into an xz axes submodel control subvector 
and a y axis submodel control vector of the form shown for 
the five body model in equations (4-5) and (4-6). 

4.4 COMPARISON OF PARTITIONED FIVE BODY AND 
SIX BODY ATTITUDE CONTROL PROBLEMS 

1 . The five body and six body xz axes and y axis submodels 
were partitioned in the manner suggested by equations (4- 
3) through (4-6) and compared as shown in the following 
equations . 
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[01 


0 


n x n state subvector coefficient submatrices for n 
body model (n = 5, 6) 

n x n null submatrices 

n x r control subvector coefficient submatrices 
n x r null submatrices 


Solid lines denote boundaries of submatrices of 6 body 
model . 

Dashed lines denote locations of rows and columns added 
due to addition of 6th body. 

From equations (4-18) and (4-19) it is evident that addi- 
tion of body 6: 

a. added four new rows and four new columns to the state 
vector coefficient matrix of the xz axes submodel; 

b. added four new rows and increased the upper bound on 
the number of columns of the control vector coeffi- 
cient matrix of the xz axes submodel by two; 

c. added two new rows and two new columns to the state 
vector coefficient matrix of the y axis submodel; 

d. added two new rows and increased the upper bound on 
the number of columns of the control vector coeffi- 
cient matrix of the y axis submodel by one. 
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2 . 


Algebraic expansions of the submatrices revealed that 
elements of the state vector coefficient matrices of the 
xz axes and y axis submodels of the five body configura- 
tion were affected by the addition of body 6 as follows. 

a. Null elements remained null elements implying that 
null submatrices remain null submatrices. 

b. Elements lying in the columns and rows associated 
with the body to which body 6 was attached experi- 
enced much larger changes in magnitude than the re- 
maining elements. 

3. From equations (4-3) and (4-4), (4-12) through (4-14) and 
equations (4-18) and (4-19) it was seen that addition of 
body 6 increased the dimensions of , the solution of 
the Riccati equation associated with the xz axes submod- 
el, from 20 x 20 to 24 x 24 and the dimensions of , the 
solution of the Riccati equation associated with the y 
axis submodel, from 10 x 10 to 12 x 12. 

4.5 REVIEW OF ALTERNATE PARTITIONINGS 


1 . The xz axes submodel was decomposed into an x axis sub- 
model and a z axis submodel to replace the solution for 
the K x Riccati matrix of dimension 24 x 24 by the sequen- 
tial solution for the Riccati matrix, K x , of dimension 12 
x 12 for the x axis submodel and the Riccati matrix, K z , 
of dimension 12 x 12 for the z axis submodel of the six 
body configuration. 

2. Additional decompositions of the xz axes submodel and the 
y axis submodel were developed analytically. 
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SECTION 5 


5.0 COMPARING ATTITUDE CONTROL EFFECTIVENESS FOR ACTUATORS ON 
TWO BODIES OF THE MODEL WITH THAT FOR ACTUATORS RESTRICTED 
TO A SINGLE BODY OF THE SAME MODEL 

5 . 1 INTRODUCTION 

The digital computer simulations of the five body approxima- 
tion of the prototype flexible spacecraft under modular atti- 
tude control reported in Chichester (5-1) and Tiffany (5-2) 
were predicated upon restriction of the control actuators to 
the central body (body 1) of the model. This restriction 
reflected the fact that control torques should not be applied 
directly to the lightweight solar wings or panels approximat- 
ed by the remaining four rigid bodies and the flexible sus- 
pension connecting them to the central body. However, the 
addition of a sixth rigid body and its associated hinge char- 
acteristics producing the six body configuration, shown in 
Figure 4-1, introduced the possibility, if not the necessity, 
of placing actuators on both bodies 1 and 6 in order to ef- 
fect attitude control of the resulting vehicle. The possible 
need for the addition of actuators on body 6 arose from the 
fact that this appendage, representing a combination of a 
second spacecraft and the deployable truss connecting it to 
the five body model of the prototype flexible spacecraft, 
involved much higher concentrations of mass and rotational 
inertias than did the two appendages representing the solar 
wings (bodies 2,4 and bodies 3,5). Two principal questions 
to be addressed in this study were the following: 

1 . Would the modular control system have enough authority to 
control the six body model's attitude effectively if the 
actuators were restricted to the central body? 
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2 . 


Would the addition of actuators to body 6 improve the ef- 
fectiveness of the modular attitude control significantly 
without generating excessive coupling effects between the 
actuators on bodies 1 and 6? 

5 . 2 PROCEDURE FOLLOWED 

1. Numerical values of the physical parameters of rigid body 
6 and three degree of freedom hinge 6 of the six body 
configuration depicted in Figure 4-1 were chosen so that 
the oscillatory modes of this single body and hinge would 
closely match the lowest oscillatory modes of the assem- 
bly of bodies and hinges representing the truss and 
spacecraft 2 in the ten body configuration shown in Fig- 
ure 2-2. 


a. The location of the center of mass and the principal 
moments of inertia about this center of mass were de- 
termined for the five serially connected bodies rep- 
resenting the truss and spacecraft 2. 

b. A two body model consisting of the single body model 
of the truss attached in cantilever fashion to a wall 
of very high mass and inertia was developed. 

c. Equations were developed for calculating the spring 
coefficients for the three degree of freedom hinge 
between the bodies of the two body model as a func- 
tion of the natural frequencies of the undamped model 
of this configuration and its rotational inertia co- 
efficients. 

d. Equations were developed for calculating the damping 
coefficient corresponding to n% damping for each 
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spring coefficient associated with the joint in the 
two body model of the truss as cantilever. 

e. The natural frequencies associated with the lowest 
oscillatory mode about each spatial axis in the five 
body model of the four body representation of the 
truss attached to a wall were equated with the natu- 
ral frequencies of the two body cantilever model. 

2. The xz axis submodel, equation (4-18), and the y axis 

submodel, equation (4-19), of the six body configuration 
were modified for restriction of actuators to body 1 and 
also to bodies 1 and 6. 


a) Restriction of the actuators to the central body 

(body 1) reduced the control subvectors, u , u and 
u y to u ix , u lZ and u iy , respectively, and the and 

0 submatrices to the dimensions 6x1. 


b) Corresponding to restriction of the actuators to bod- 
ies 1 and 6 the control subvectors were truncated to: 


U = (U , u ) ; 

—X lX (, x 


u = ( u , u ) ; 

-y iy' &y 


H z - <v- V> • 


The second through fifth columns were then removed 
from the B^j and 0 submatrices reducing their dimens- 
ions to 6 x 2. 


3. The xz axes submodel and the y axis submodel for the six 
body configuration were simulated on the digital computer 
for restriction of actuators to body 1 and for restric- 
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tion of actuators to bodies 1 and 6 under modular atti- 
tude control. 

4. Initial angular displacements of .5 degrees with respect 
to each spatial axis were applied to the central body 
(body 1) and at each of the three degree of freedom hing- 
es between the rigid bodies of the model and the result- 
ing responses were plotted and tabulated. 

5.3 RESULTS 


The maximum magnitudes of the Euler angles of the rigid bod- 
ies comprising the xz axis submodel with selected distribu- 
tions of actuators and sensors are listed in Tables 5-1 and 
5-2. Since the magnitude of each initial angular displace- 
ment was .5 degrees, entries in the tables of substantially 
greater magnitude imply considerable overshoot in the re- 
sponses. Each row of either table represents the maximum 
magnitudes of the Euler angles of a spatial axis for the set 
of conditions listed at the left side. The angles with sub- 
script "1" are the inertially referenced attitude angles of 
the central body (body 1). The double subscripts identify 
the interface with which the corresponding Euler angle is 
associated (4> 1Z refers to the angle about the x axis between 
bodies 1 and 2). The angle labels at the top of the table 
are arranged in approximate spatial order as depicted in 
Figure 4-1. 

From these tables the following were observed. 

1 . The largest maximum Euler angle magnitudes for a given 
set of conditions generally occur at the outer extremi- 
ties of the solar panels. 
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TABLE 5-1 

X AXIS MAXIMUM EULER ANGLE RESPONSES FOR 
DIFFERENT DISTRIBUTIONS OF ACTUATORS AND SENSORS 


U1 

i 


Actuators and 
Position Sensors 
on Bodies 1 & 6 

^24 & ^35 Sensors 

Actuators and 
Sensors on Body 1 

^24 & ^ 35 Sensors 
with Body 6 Sensors 



35 1 max 

1 ^ 13 1 max 

llmax 

[deg.] 

[deg.] 

[deg. ] 

1.16 

.50 

.612 


1.18 *50 .627 


1.24 .50 .50 


16 [max 

l^^lmax 

1^24 1 max 

[deg.] 

[deg.] 

[deg.] 

.50 

.50 

1.12 


.758 .50 1.13 


.514 .50 1.55 



1 1 1 ] 1 ] ] 1 1 


cn 

1 

ON 


TABLE 5-2 

Z AXIS MAXIMUM EULER ANGLE RESPONSES FOR 
DIFFERENT DISTRIBUTIONS OF ACTUATORS AND SENSORS 

|*35|max | v 13 Imax |*llmax l*16lmax l*12lmax 

[deg.] [deg.] [deg.] [deg.] [deg.] 


Actuators and .849 .50 .688 

Position Sensors 
on Bodies 1 & 6 

^24 & *35 Sensors 

Actuators and Position 
Sensors on Body 1 

^24 & ^35 Sensors 

with Body 6 Sensors .856 .50 .725 


.572 .50 


.732 .50 


Actuators and 
Position Sensors 
on Bodies 1 & 6 


*24 & *35 
Sensors with 


1.10 .50 


.507 


.50 .50 


1*24 Imax 
[deg.] 

.849 


.856 

1.12 


Asymmetric 
Disp lacement s 



2. Elimination of the control actuators from body 6, one of 
the two bodies on which they can be placed, increases the 
overshoot in the responses of the central body (body 1) 
and the responses of the relative Euler angles between 
bodies 1 and 6 to the point where they are very poorly 
controlled. 

3. Application of asymmetric initial displacements in the 
Euler angles at the three degree of freedom hinges be- 
tween the rigid bodies of the model produced Euler angle 
responses with especially high overshoots close to the 
outer extremities of the solar panels. 

Responses of the relative Euler angle between bodies 1 and 6 
for the y axis submodel with selected distributions of ac- 
tuators and sensors are presented in Figures 5-1 and 5-2. 
Since this submodel is linear, the responses actually are 
sinusoidal in nature although, in some instantances, the 
characteristics of the plotting subroutine utilized may 
result in their appearing to be comprised of contiguous 
straight line segments. 

From Figures 5-1 and 5-2, the following were observed. 

1. Elimination of control actuators from one of the two bod- 
ies on which they could be placed, body 6, introduces an 
offset in the final or steady state value of the respons- 
es of the relative Euler angle between bodies 1 and 6. 

2. Elimination of control actuators from body 6 renders con- 
trol of the relative Euler angle between bodies 1 and 6 
completely ineffectual about the y axis even with sensors 
present on body 6 which was seen by comparing Figure 5-2 
with Figure 5-1 . 
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SECTION 6 


6.0 CONCLUSIONS AND RECOMMENDATIONS 


During the period covered by this report the digital computer 
simulation of the ten body approximation of the three axis 
rotational dynamics model of a prototype flexible space plat- 
form consisting of two spacecraft interconnected by the 
MSFC/hybrid deployable truss with modular control was accom- 
plished except for the completion of validation runs. An 
important part of this task was the generation of numerical 
values of the physical parameters of a five body approxima- 
tion of the three axis rotational dynamics of the MSFC/hybrid 
truss in cantilever configuration that matched five modes of 
the dynamic responses of this approximate model with those of 
a finite element NASTRAN model of the truss that was sup- 
pi ied . 

The remainder of the effort during this time period utilized 
a five body approximation of a three axes model of the rota- 
tional dynamics of a prototype flexible spacecraft and a six 
body approximation of the ten body model of the flexible 
space platform described above. The five body model was 
utilized to evaluate the effects of elimination of state 
variable sensors upon the modular attitude control of the 
prototype flexible spacecraft and to generate sensitivity 
coefficients of its state variables with respect to the mod- 
el's parameters. The six body model was used to evaluate 
software changes required to accommodate the addition of 
another body to the five body model and to compare the effec- 
tiveness of attitude control for actuators on two of the 
model's bodies with that when the actuators were restricted 
to a single body of the same model. 
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6.1 CONCLUSIONS 


The following conclusions are based mainly upon modeling and 
digital computer simulation of a four body approximation of 
the three spatial axes rotational dynamics of the MSFC/hybrid 
deployable truss, modeling simulation and control of a five 
body approximation of a prototype flexible spacecraft and 
modeling, simulation and control of a six body approximation 
of the space platform which consists of two spacecraft inter- 
connected by the MSFC/hybrid deployable truss. 

1 . Five modes of the rotational dynamics of a three dimen- 
sional finite element NASTRAN model of the MSFC/hybrid 
deployable truss were closely approximated by matching 
the eigenvalues of a four body model of the truss with 
those of the NASTRAN model. 

2. The close approximation of the dynamic modes of the 

NASTRAN model of the truss by those of the four body mod- 

el supported the incorporation of the latter model in the 
ten body model representing two spacecraft interconnected 
by the MSFC/hybrid deployable truss. 

3. Due to the light damping present in the vehicle model, it 

was difficult to control oscillations at the free ends of 

the appendages with control actuators restricted to the 
central body. 

4. Increasing the magnitude of the elements in the diagonal 
weighting matrix, W u , for the control vector in the quad- 
ratic performance index tended to decrease the amount of 
control torque expended about a given axis but it also 
reduced the effectiveness of attitude control about that 
axis. 
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5 . 


Increasing the magnitude of the elements in the diagonal 
weighting matrix, Q, for the state vector in the quadrat- 
ic performance index tended to decrease the peak magni- 
tude of the corresponding state variable but it also re- 
duced the effectiveness of the attitude control about the 
corresponding axis. 

6. The changes in parameter values of the five body model to 
which the scalar state variables (axial components of 
body angular rates and Euler attitude angles of the rigid 
bodies in the model) are most sensitive, are in descend- 
ing order, changes in the z axis damping coefficients, 

c s bb z' c s^z' C S13Z' c sizz' and changes in the z axis 
components of the location vectors of the hinges connect- 
ing the central body (body 1) to the adjacent bodies, 
v l l z r and V 1ZZ- 

7. Among the seven scalar state variables of the five body 
model showing the highest peak sensitivity and the seven 
state variables showing the highest steady state sensi- 
tivity, six parameters are common to both groups. 

8. The computer software developed for the five body model 
may be extended to accommodate the addition of another 
body if provisions are made for the expansions in the 
dimensions of the subvectors and submatrices of the par- 
titioned form of the model described in equations (5-18) 
and (5-19) and for the changes in the values of the ele- 
ments of the matrices, especially those associated with 
the rows and columns corresponding to the rigid body to 
which the sixth body is attached. These concepts also 
may be applied to the case in which another body is added 
to the ten body model of the space platform. 
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9. For the six body approximation of a flexible space plat- 
form, depicted in Figure 5-1: 

a. Control actuators and sensors must be present on both 
bodies 1 and 6 and sensors of relative Euler angular 
displacements about the z axis close to the extremi- 
ties of the solar panels are required for effective 
three axis attitude control of the spacecraft in re- 
sponse to symmetric initial angular displacements. 

b. Elimination of the control actuators from body 6 left 
insufficient control authority for effective attitude 
control even with the appropriate sensors on body 6. 

6.2 RECOMMENDATIONS 

The following items are recommended for future study concern- 
ing control of flexible spacecraft. 

1. The effects of incomplete state feedback upon the atti- 
tude control of the five body model of the prototype 
flexible spacecraft should be investigated. 

2. The partitioning applied to the five body state variable 
model of the flexible spacecraft should be extended to 
the Riccati equations associated with the corresponding 
optimal attitude control problem. 

3. A prototype sigle axis model of a flexible spacecraft 
should be developed to facilitiate comparison of differ- 
ent attitude control approaches with minimal computer 
requirements . 
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4. Procedures should be developed for estimating the real 
time computational requirements for the controlled space- 
craft models. 

5. The attitude responses of the flexible spacecraft models 
with modular control should be compared with those ob- 
tained for the same models with classical frequency do- 
main control to determine those conditions under which 
modular control might be advantageous. 

6. A control concept consisting of determination of the num- 
ber of control moment gyros, determination of a control 
law and design of a mission operations momentum manage- 
ment plan should be developed for at least one represen- 
tative spacecraft. 

7. Methods of momentum management for counteracting the bias 
portion of the external torques on the space vehicle 
should be developed. 

8. A momentum management algorithm should be developed under 
which a spacecraft automatically adjusts its attitude 
with respect to the orbit plane in such a way as to mini- 
mize the angular momentum buildup. 
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